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Abstract 

The photon self energy tensor in a crystal with double exchange and ferromagnetic 
or canting magnetic orders is analyzed at zero temperature with emphasis on the signal 
of the time-reversal invariance breaking. Three comparable contributions, spin suscep- 
tibility, local magnetic field and spin-orbit coupling arc estimated and the prospect of 
their experimental detection is discussed. 



1. Introduction 



Recently, there are growing interests in the double exchange mechanism 
proposed by Zener [1], Anderson and Hasegawa [2] over forty years ago. It 
stems from the possible relevance of the mechanism to the colossal magneto- 
resistance of the single perovskite material La2- x Ca x MnOs [3] and is con- 
jectured also to be responsible to the newly synthesised double perovskite 
material (Sr, Ba)2Y Ru\- x Cu x Oq [4]. The physical picture of the double 
exchange mechanism in the latter system lies in the outer electronic con- 
figuration of the -Rit-ion. Prior to doping, the Ru-ion carries a charge +5 
and its outer electronic configuration is 4<i 3 , the spin of each electron are 
forced parallel by strong Hund's rule coupling and the three electrons fills 
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in a closed multiplet of the cubic group. The total spin of each Ru +5 is 3/2 
and interacts with each other through an antiferromagnetic super-exchange. 
The parent compound is therefore an antiferromagnetic insulator. When the 
charge balance is offset by doping of lower valence transition elements, Cu, 
holes are introduced to the outer electronic configuration and they start to 
hop. The system becomes a conductor and may even be a superconductor 
below certain temperature. The spins of these holes couple strongly with the 
Ru +5 ions via Hund's rule and bring in ferromagnetic couplings between ion 
spins. The prototype Hamiltonian of this mechanism reads 



H — ^ Ji jSi ■ Sj \i c ^ Bj ■ Sj 

ij 3 

~\ + *J*<) + ^ A E + Hcoui.. (1.1) 

ij 3 

where Sj is the total spin operator of the ion core, ty - are the creation and 
the annihilation operators of the itinerant carriers, each of two components 
with spin up and spin down, Jy denotes the antiferromagnetic exchange, 
Uj denotes the hopping amplitude of the itinerant electrons and A denotes 
the Hund's rule coupling of the itinerant electrons with the ion cores. The 
Coulomb interaction, Hcoul.i is given by 

„2 , ■ vl>! vCr- vCr^ ,vl> . , 

(1-2) 

with the first term the on-site Coulomb repulsion and e the dielectric constant 
produced by the ion cores. 

The competition between the antiferromagnetic exchange and the fer- 
romagnetic coupling may lead to a canting magnetic order below certain 
temperature, as was analyzed by de-Gennes [5] and supported by specific 
heat measurement of the double perovskite materials [6]. In the previous 
paper [7] , hereby referred to as paper I, we studied the magnetic ordering of 
the Hamiltonian (1.1). Spinwave spectrum was calculated from the Hamil- 
tonian (1.1) with a canting magnetic order. The important role played by 
long range Coulomb interaction in maintaining such an order is emphasized. 
We also performed extensive Monte Carlo simulations on de-Gennes's effec- 
tive magnetic Hamiltonian and computed the specific heat as a function of 
temperature. All our results were found consistent with the experimental 
measurements of the specific heat and the magnetic susceptibility. 

In this paper, we shall address to the optical properties of the unusual 
magnetic ordering caused by the double exchange. The characteristic term 
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of the Hamiltonian (1.1) is the Hund's rule coupling term. A reasonable 
estimate [2] gives rise to XS ~ leV. In the paramagnetic phase, such a term 
causes strong scattering of itinerant holes. In the ordered phase with ferro- 
magnetic order or canting, such a term represents an effective homogeneous 
magnetic field acting on the hole spins only. With the magnetic moment of 
a free electron, 2 //Bohr; the equivalent magnetic field strength is as high as 

S equiv ~ 10 4 Tesla, (1.3) 

(2 / UBohr-^equiv ~ XS). Detectable effect signaling the violation of time reversal 
invariance is expected and such effects are the main theme of this work. 

In the following two sections, we shall couple the system described by the 
model Hamiltonian (1.1) with radiation field calculate the photon self energy 
function of the medium. A f.c.c cubic structure like the double perovskite 
materials and and a tetragonal layer type magnetic ordering are assumed for 
concreteness. In contrast to the ordinary magnetooptical effect, the equiva- 
lent magnetic field (1.3), essentially rooted in the Coulomb forces, does not 
acts on the orbital motion directly. The optical activity comes from three 
different sources: 1) the Pauli term; 2) the secondary magnetooptical effect 
because of the real magnetic field produced by ordered spins; 3) a spin-orbital 
coupling which comes together with the Pauli-term in the nonrelativistic ap- 
proximation of the Dirac electrons. The prospect of the experimental de- 
tection, the relevance of possible genuine superconductivity order and the 
validity of various approximation we made are discussed in the final section. 



2. The Electromagnetic Coupling and Its Parabolic Approximation 

The Hamiltonian combining the antiferromagnetic exchange, the double 
exchange and the electromagnetic coupling reads 

H=f d 3 r±(El + B 2 )+Y1 JiA ■ Sj 4" • $j 

ij j 

+ 1 5^ ^ ^ + ^ t^7i ) - ^ A ^ 5^ - M ^ ^ ^ +- ^ - ^ ^ ^ -h ^Oo^Z . , 

ij j j 

(2.1) 

The Coulomb gauge is adapted with V • A = 0. The transverse electric field 

Etr. = — Off and the magnetic field B = V x A where Bj = B(Rj) with Rj 
the position vector of the jth site. The lattice gauge coupling to the carrier 
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field is given by 



Uij = exp 



ze 



[ dslij • A(Ri + slij)] (2.2) 
Jo J 



where kj denotes the displacement vector from the ith site to the jth site. 
In what follows, we shall restrict the hopping amplitude within the nearest 
bond only. 

Following the conventions in I, we introduce the lattice frame (x, y, z) 
with x, y and z parallel to the (100), (010) and (001) directions of the crystal 
respectively. The locations of magnetic ions (total number = Af) are given 
by 

R = n\e\ + niei + uses, (2.3) 
where ni, ni and 773 are integers, 

ei = ^(y + z), (2.4) 

e 2 = l -{z + x), (2.5) 

e 3 = ^(x + y), (2.6) 

with / the side length of a basic cube (/ = 8.16A, for Sr2Y ' Ru\^ x Cu x Oq). 
Adapting the tetragonal layer magnetic ordering, we divide the lattice into 
two sublattices shown in Fig 1. with 77,1 + 77-2 = even on the sublattice A and 
n\ + 77-2 = odd on the sublattice B. The expectation value of the ion spins 
of the sublattice A is denoted by Sa = S(a and that of the sublattice B by 
S b = S(b- We choose the £-axis along the direction of Sa + Sb and the 
77-axis as 

« = ^ (2-7) 
sin r) 

with © the mutual angle between Sa and Sb- We introduce further 

i = V x C, (2.8) 

The set (£, fj, Q forms a right hand coordinate system and will be referred to 
as the spin frame. Without the magnetic anisotropic energy, the orientation 
of the spin frames does not couple with the lattice frame. 

With large spin approximation, the Hamiltonian (1.1) can be written as 

H = H em _ + H sw _ + #Coul.> (2-9) 
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(a) (b) 

Figure 1 (a) The magnetic ordering. The magnetic sites of the sublattices A and B are 
marked explicitly, (b) The spin frame. The angle between (a, and (b axes is G. 

where H e m is the part of the Hamiltonian we are concerned with. Its is 
given by 

H . m . = j d 3 f^(El + B 2 ) -\tJ2 (<U ab ^ b + *lu ba y a ) 

J ~ <ab> 

E (^aU aa ,y a , + yl,U a , a y a )+ ^(*lu w * v + 9\,U Vb * b j 

<aa'> <bb'> 
a b a b 

(2-10) 

where t and t' denote the interlayer and intralayer hoppings respectively. The 
terms containing spin wave operators are included in H s w and have been 
dealt with in I. 

Within Heisenberg representation, the time reversal operator T trans- 
forms 

TA(r, t)T- 1 = -A(r, -t), (2.11) 
T^jifyT- 1 = K(j 2 ^j{-t) (2.12) 

and 

TtfJ^T" 1 = K*v](-t)a 2 (2.13) 
with k a phase factor. It is easy to see that 



TH e . m T~^ H e . m ., (2.14) 
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because of the terms containing Pauli matrices and therefore -£Ze.m. is not 
invariant under a time reversal. 

We may also introduce a pseudo time reversal operator T, whose repre- 
sentation is 

TA(r, t)T- x = -A(f, -t), (2.15) 
TV^T- 1 = n%(-t) (2.16) 

and 

T^t^T" 1 = «/*#](-*). (2.17) 

with k' another phase factor. It is interesting to notice that except the last 
two terms of (2.9), the rest of -fie.m. is invariant under T. This property is 
very useful for the one-loop calculations in the next section. The Hamiltonian 
(2.1) is also invariant under the space inversion 

PA(r, t)P ] = -A(-r, t) (2.18) 

and 

P^j{t)P ] = r}Vj*(t) (2.19) 

with j* the inversion image of j and rj a phase factor. 

The complicated dependence of U\j makes the perturbative expansion in 
terms of the fine structure constant rather cumbersome, we shall make a 
parabolic approximation for the rest of the paper. We found in I that, at 
the absence of the vector potential, (2.10) reduces to J^p^pEp^p w ^ n 4x1 



column matrix typ and Ep the 4x4 hermitian matrix. 



where 



and 



with 



Ep= -t'up-tvppi - -\S(p+( A ■ t + p-( B ■ r), (2.20) 



n = ( n' °U=(? "Wl? 1) (2.21) 



(7i I ' A \ a 2 r 6 \ £7 3 



\p a ,n] = 0- (2-23) 
The kinetic energy, ^pty^Epfyp, is invariant under the C-transformation, 

C^pC^ 1 = —p lTs^p and C 2 = 1. This symmetry can be diagonalized through 
the unitary transformation typ = Vi^p with 

V = -j= (p- + P+t 3 ) (pi + p 3 ) ( cos j - ir 2 sin ® ) . (2.24) 
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We find that 

V^EpV = -t'up - tv^ps (V 3 cos ^ - n sin ^ - ^XSr 3 , (2.25) 

which is easily diagonalized and produces four energy bands: 

e p =-t'(up+n)±A± (2.26) 

with 

A± = ^\ 2 S 2 ±\Stvpcos®+t 2 v 2 , (2.27) 

where ± of Ap labels the eigenvalue of the C operator. The lowest band 
corresponds to 

ep= -t'iup+n) - At (2.28) 

which overlaps partially at the top with the next lowest band, i.e., ep = 
—t'(up+p,) — AZ. In what follows we shall assume that only the band (2.28) is 
partially filled, a likely situation with doping conditions of (Sr, Bd)2Y Ru\- x Cu x Oq^ 
With the parabolic approximation we have 

up=2-j(p 2 x +p 2 y ) (2.29) 

and 

f 

Vp = 4--(p 2 x +p 2 y + 2p 2 z ). (2.30) 

Regarding iftp as the Fourier component of a 4 x 1 fermion field ip(r) in 
continuum, we have jj J2p = I (f^s an d the " coordinate" representation of 
kinetic energy operator, (2.25), reads 

/l 11 1 (9 $ 

L2 \m m/ moz oz 

-v^pxr'^ - pc^tp - -XSip^rsip , (2.31) 

where the parameters m, m' and v are related to the discrete Hamiltonian 
through m'- 1 = t'l 2 /2, mT 1 = tl 2 /2, v = 4t, \x c = A* + 2 *'- Replace the 
ordinary derivatives by the covariant ones, we obtain the continuum approx- 
imation of (2.10), i.e. 

H e . m . = [ d 3 fhE 2 T +B 2 )+ [ d3f[^(V ± +iel L )V' t -(^7+-P3^)(V ± -iel L )J 
J 2 J 12 \m m / | 



+— + ieA z )^ p 3 T^(-^- - ieA z )ip 
m oz oz 

-u^psT^tfj - iJL c ^i> - -\S^ni> - y^s^tt • Bijj , (2.32) 

where 

a = V^tV = -pn{i - piT 2 f] + r^C, (2.33) 

7-3 = r 3 cos ® - rising, (2.34) 

. e e . 

t x = 73 sin — + Ticos— , (2.35) 

e is positive for holes and negative for electrons and we have substituted 
g e u B for fj, h . 

It is well known that, the Hamiltonian (2.10) and (2.32) produce the 
same physics as long as the length scales involved are much larger than the 
lattice spacing, /. The Hamiltonian (2.32) is not bounded from below but 
this does no harm for the perturbative calculations provide the condition 
of the parabolic approximation is met. Like the Hamiltonian (2.10) the 
Hamiltonian (2.32) is not invariant under a time reversal. But the pseudo 
time reversal transformation leaves all terms except the last one (Pauli term) 
invariant. 

The electric current operator of the system consists of three terms 

J(f, i) = J d (r, t) + J p (r, i) + J s (f, t) (2.36) 

— > — » — » 

where J p and J s denote the diamagnetic, paramagnetic and spin part 
and their expressions are listed below 

2 2 

Jd(r,t) = ~^(— + -A*t$W± - -^P3T^A z z; (2.37) 
2 \m m J m 

Mr, t) = -f [*t (-L + I^) - h.c] - I (*W 8 ^ - h.c.)i 

(2.38) 

and 

J s (f,f) = |^Vx^. (2.39) 
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Figure 2 The Feynman diagrams for the photon self energy tensor. 
3. The Response Functions 

When the system is illuminated by a beam of light with a wave vector k. 
The optical response is entirely determined by the transverse components of 
the retarded photon self energy tensor in the light frame, i.e. 

a a p{u) = (e a )i n §(w, k)(e a )j (3.1) 
where, at zero temperature, 

ReU^(oj, k) = ReUij(u, k) (3.2) 

and 

Imng(w, k) = sign(w)Imn ij (w, k) (3.3) 

with II(u;, k) time-ordered self energy tensor which can be calculated dia- 
grammatically. For a system invariant under the space inversion, we expect 

U ij (oj,k) = U ij (oj,-k). (3.4) 

For a system invariant under time reversal, we have 

Uij(u, k) = Uji(u!, -k). (3.5) 

Come to the system we are considering, the time reversal invariance is spon- 
taneously broken by the magnetic ordering but inversion symmetry is intact. 
The self energy tensor will acquire an antisymmetric part, i.e. 

Uij(uj,k) -Uji{u,k) ^0, (3.6) 

which we shall focus our attention on for the rest of the article. 

The one-loop diagrams for Uij(u, k) are displayed in Fig. 2. 

The evaluation of them is straightforward and we tabulate below the 
results under the approximation XS » t,t' u> ~ XS and p/j/ 3 << 1 where 
Ph is the number of itinerant holes per unit volume, 

/(ftp 
j^s°(Pc - ep) (3.7) 
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with 9(x) = 5(1 + sign(x)). We divide Hij(u> : k) into three terms, 



Ily (a;, k) = U\f(u, k) + U].'(^ k) + njJ'( W) k) 
The IlW term comes from the orbital motion only and is given by 
u\f(oj, k) = e 2 p h —,{8ij - 5 iz 5 jz ) + — cos ^(Sij + 5 iz 5 jz ) 



r(2), 



r(3), 



.m 



m 



-i f d 3 fe~ ilr f dte iu)t < \TJ p j(r, t)J p j(0, 0)| > 
J J — 00 



e 2 p fe AS 
m 2 A 2 **? 2 — uj C/L 



- SizSjz) + — COS ^r{Sij + 5 iz 5 jz ) 



m 



m 



(3.8) 



<P±> - &iz&jz) + 8 <p 2 z > 5 iz 5 jz sin 2 (3.9) 



where 



The n^ 2 ) term of (3.8) represent the cross term between the paramagnetic 
current and the spin current and is given by 



U^\u,k) = -i / cfre- lt? 



dte wt < \T[J Pji {r,t)J sJ {0,0)+J Sji {f,t)J pJ (0,0)}\ >\ 



'i[ e imnkmCn^-j ( w ? k) ej mn k m (^ n Ai(u , k)] 



(3.11) 



where 



Aj(u, k) = g e 



VBeph 
2mu 



( m © 
— 7 cos — 
\rri 2 



6 A 2 S 2 cos 2 f -cA . x . 
+ A 2 S 2 -. 2 



A 2 S 2 cos 2 f - W 2 



A 2 S 2 -^ 2 »z~3*y ( 3 - 12 ) 

Finally, the II^ 3 ) term is related to the spin susceptibility tensor Xij{ u ) y i a 



u\f(u,k) = / dte luJt < \TJ 8i i(f,t)Jsj(P,0)\ >= t imn e jm i n ik m k m iXnn'(^), 
J — 00 

(3.13) 

Xij(v) = Riv R jvX»v(u) (3.14) 
10 



where 



with the indices n, v referring to axes of the spin frame (p, v = £, rj, Q and 
Rip the elements of the orthogonal matrix 

(x • £ x • r) x ■ O 
y-{ y-v y-C • (3.15) 

The nonvanishing components of x Q( a(w) are 

1 2 2 AS cos 2 ! 



a,' COS y 



*CCM = ~^eP^I A2| g2 _^2 Sin2 f • ( 3 - 19 ) 

We observe that 1) The antisymmetric part comes only from II^ 2 ) and 
terms (the statement is true without the approximation AS 1 >> t, £' and 
n/j/ 3 << 1), while remains symmetric, a consequence of T-invariance. 
2). With the large double exchange approximation, the symmetric part of 
riy(u;, k) is dominated by the first term of (3.9) and the antisymmetric terms 
are smaller than this term by a factor of the order of k 2 /m\S. For to not at 
the pole, ±XS, we have Ilfj(uj, k) = Hij(u, k) and from this the expression 

of a a p{u, k) follows. 

The propagation mode of a photon is determined by 

[(u 2 - k 2 )S aP - a af3 (u, k)\<j>p = (3.20) 

with (j) the polarization vector, ± k. The T-violation implies that the 
eigenmodes are elliptically polarized instead of linearly polarized and the 
system possesses the optical activity. 

Consider a special circumstance, in which k, ( and z are all parallel, 
and the magnetic order is purely ferromagnetic. The self energy tensor <j a p 
reduces to 

<r af} (u, k) = e 2 Ph (^j + ~y a p - ^ 2 e Q/3 ff e W| A 2^_^2 ' ( 3 - 21 ) 
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where we have shown explicitly the symmetric and antisymmetric parts of 
a a fj(uj, k) and only the leading order terms are retained in each part. The 
eigenmode equation (3.20) are reduced to 

(iV 2 - UJ 2 )5 a p - n 2 UJ 2 (l - ^apglphl^B X 2 S 2 _ ^2 ) = °> ( 3 - 22 ) 

where we have introduced the index of refraction, n(u) via k = n(uj)u and the 
plasma frequency u 2 = e2 Ph(^/+^ cos §)• Clearly, the eigenmodes of (3.20) 
are two circularly polarized photon with opposite senses. The corresponding 
indices of refraction differ by an amount 

An(w) = \V^9Pb X 2 S 2_ U 2 ( 3 - 23 ) 

with eo = 1 — which gives rise to a Farady angle 

A6 = nAn{u) (3.24) 

per wavelength. For Sr2Y Rui- x Cu x Oq, we have ph = 2x/v with v = Z 3 /4 = 
136A 3 . Assuming g e = 2, = and ^- + \ = -J-, we obtain: 

lib lib I lb(^ 

A6 = 0.62 x lO~ A ,y^ XUJ 9 (3.25) 

with XS, ui in ev and AO in radians. 

There are two additional contribution to the T-violation signal which are 
not included in the Hamiltonian (2.32) and their magnitudes turn out to be 
comparable to (3.25). 

1. The secondary magnetooptical effect The aligned ion-spins will pro- 
duced a net magnetic field, which gives rise to the ordinary magnetooptical 
effect [8]. For a magnetic field B = Bz and k = kz, the Farady angle reads 

AO' = tt-^ (3.26) 

with ub = e B{^r + ^ cos ®) the cyclotron frequency (See appendix A for 

the derivation with a general anisotropy). For a spheroidal sample with M 
parallel to its symmetry axis, the average magnetic field is 

B = 7M 

with M the magnetization given by 

M = cos — 

v 2 
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with gi the Lande factor of the core ion, where the contribution from the 
itinerant holes has been neglected on account of low dopping condition. The 
factor 7 changes from zero in the limit of a thin disc perpendicular to the 
field to one in the limit of a long rod along the field. Assuming gi = 2, S = §, 
= and + ^ ~ we find the corresponding Farady angle 



A6' = 1.89 x Kr 3 ^-^. (3.27) 



e ur 

Though this effect could be considerably larger that (3.25) but the double 
exchange effect may still be detectable because of its resonance nature. 

2. Spin-orbit coupling It is interesting to know there is a relativistic term 
which turns out to to comparable with (3.25). The physical reason is the 
same as that for the spin-orbital coupling of atomic levels. In what follows, we 
give a classical argument [9] and defer the quantum mechanical treatment to 
the appendix B. For an electron(hole) moving in an external electromagnetic 
field (E, B) at velocity v(v << 1), The spin precession equation reads 

^ = sx (— B'-ujt), (3.28) 
dt m e ' y ' 

where B' is the magnetic field in the rest frame of the electron(hole) and oot 
is the angular velocity of Thomas precession. We have set g e = 2 in (3.28) 
since we do not know how to renormalize the bare value in a medium for this 
effect. The equation (3.28) can be derived with a Hamiltonian term 

H' = -(2e^ B B' - Cut) • s. (3.29) 
To the first power of the velocity, we have 

B' = B — v x E (3.30) 

and 

Q T = --vx a = -/j b vx E (3.31) 

with a = -^qE to the order of the approximation. Substituting (3.30) and 
(3.31) into (3.29) and express v in terms of the canonical momentum and 
the gauge potential, v = ^(p — eA), we obtain 

H' = -2^ B s- B - —8- E x (p-el) (3.32) 
2m P 
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as an improved Pauli term. The term E x A pick up the information of the 
violation of the time reversal invariance with the corresponding Farady angle 
per wavelength given by 



A6" = 



™ 2 pI e 



cos — = 0.62 x 10 



4 X 







(3.33) 



4. Discussions 

In the previous sections, we analyzed the photon self energy tensor in 
a canting magnetic ordering, focusing the signal of the T-violation. There 
are three contributions to the optical activity, eqs. (3.25), (3.27) and (3.28). 
Although the equivalent magnetic field due to the double exchange is huge, 
the magnitude of the Farady angle turns out rather small. The reason for this 
is that the double exchange acts only on the carrier spins. The orbital motion, 
neglecting the weak magnetic field produced by the ordered ion spins, does 
not give rise any T-violation signal, on account of the T invariance defined 
in eqs. (2. 15-17). Such an small effect lies at the border of the detection 
precision for a sample of one wavelength thick, as inferred from the delicate 
circular dichroism experiment, designed to test the anyon model of the high 
Tq superconductivity [10]. On the other hand, the Farady angle associated 
with T-violation accumulates as the photon traverses back and forth inside 
the sample [9]. The effect may be enhanced through multiple reflections at 
the inner surface of the sample. 

For the double exchange energy XS comparable with the photon energy, 
u, resonance character will show up in the contribution from the spin sus- 
ceptibility, (3.25). A closer look will reveal an interesting point which has 
been watered down in the leading order of the strong double exchange ap- 
proximation made in the eqs. (3.9-19). According to the Hamiltonian (2.32), 
the transition involved in the photon self energy includes both C preserving 
ones and C flipping ones. In the former case, the energy denominator for the 
transition reads 




(4.1) 



while in the latter case, the corresponding energy denominator is 




(4.2) 
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Upon integration over p within the Fermi sea, the resonance signal will be 
slightly smeared, such a smearing effect for the C flipping transition is of one 
order higher than that for the C preserving transition. The transition con- 
tributing to x^(^) is entirely C flipping type and the resonance signal in the 
Farady angle (3.25) should be fairly sharp. This feature is also independent 
of the parabolic approximation we made in the calculation. 

The frequency dependences of the other two contributions are rather 
robust. The former involves the transition among different Landau levels 
caused by the real local magnetic field of the ordered magnetic moments of 
the ion cores. This field is in general fairly weak, of the order of 1000 Gauss, 
the level spacing associated with this magnetic field is tiny in comparison 
to the photon energy, u> and the energy denominator is dominated by the 
photon energy. The relativistic term (3.28), however, is the consequence 
of the virtual transition between the Dirac sea and the Fermi sea and the 
corresponding energy denominator is dominated by twice of the rest energy 
of an electron. 

Experimentally, superconductivity was also discovered the double per- 
ovskite samples. A natural question arise as to the effect of the long range 
order on the optical activity calculated with normal electrons. For the pho- 
ton in the visible region, its energy are much higher than the gap-energy 

(2) ~* (3) ~* 

and correction to and XIV - ; (o;, k), (u, k) is of higher orders. Only thing 

we need to watch out is whether there is small term coming from U^\uj, k). 
The long range order can be introduced by adding a pairing term to the 
Haniltonian (2.32). The total Hamiltonian reads now, 

<? 

where ap, at stands for the annihilation and creation operators of the itiner- 
ant holes, and fig, fit stands for composite boson operators which represent 
Cooper pairs, i.e., 




with Yl'p extends over half of the BriUouin zone only. The pairing wave 
function, gp^, is odd under the space inversion, i.e. 

9p,q = -9~p,q (4-5) 

on account of the anticommutation relation among ay's. The superconduc- 
tivity is implemented through the condensation of the pairing operator of 
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zero total momentum, i.e. 

< S\p ?=0 \S >= VJfB, (4.6) 

where B is the long range order parameter. The wave function of the ground 
state is then 

\S >= \[{cos9p- e-^sin^V p .)|0 >, (4.7) 
v 

where cos 29 p = ^ and sin 29 p = ^ with Sp = ^/e?.+ \5p\ 2 and |%|e* 7 P = 

Bgpp -1^=0- Evidently, both the standard time reversal T and the pseudo 
time reversal T do not leave the ground state (4.7) unchanged. On the other 
hand, The following antiunitary transformation 

T'apT'- 1 = ie^ap (4.8) 

leaves both the pairing Hamiltonian as well as the ground state (4.7) invari- 
ant. This, together with the invariance under the modified space inversion in 
I, V = Pe% Q , implies that u\f(uj, 0) = U^\(j, 0) and the T- violation signal 

is suppressed relative to 11^ (a;, 0) at least by a factor k 2 /p F with p F the 
Fermi-momentum. The very fact that they must vanish as the gap energy 
vanishes pushes it much smaller than the magnitude of Xf»j(u>)- Therefore the 
estimate of the T-violation signals (3.25) and (3.28) with normal electrons 
applies to the superconducting phase as well. 

As was pointed out in [11] in the context of Raman scattering, the way 
we introduce the gauge invariant coupling, eqs. (2.11) and (3.32), requires 
that the photon frequency in the energy denominators for the transition to 
other bands not covered by the model Hamiltonian (1.1) can be dropped. 
This may be rather marginal for real system. Also the impurity effect, which 
may be relevant to the real system has been neglected. 
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Appendix A 



For the ordinary magneto-optic effect, we consider only the conduction 
band. The effective Hamiltonian with the parabolic approximation reads 

where we have introduced general anisotropy and the symmetry of (A.l) is 
orthohmobic. In our case, which is tetragonal 

1111© 
— = — = — + - cos- (A.2) 
mi 1712 m m 2 

and 

1 2 /A o \ 

— = _ cos _. (A.3) 
m^ m 2 

The gauge potential A(r, t) is chosen to produce a static and homogeneous 
magnetic field, 

B = Vx A(r, t) (A.2) 
and a time dependent and homogeneous electric field, 

The electric current operator reads 

Ja = V T (« ieA a )iP + h.c. (AA) 

m a ox a 

with a = 1, 2, 3, and the charge density operator is 

J = e^ip. (A5) 
Let | > denote the ground state of (A.l). It follows from the translation 

— * — * 

invariance of E and B and the gauge invariance of (A.l) that the expec- 
tation values < \J(r,t)\ > and < \Jo(r, t)\ > are coordinate independent. 
Introducing the operators 

1 r A 



J= - I d"fJ(f, t) (A6) 
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and 



we have 

j=<\J(r 1 t)\>=<\J\> (A8) 

and 

ep h =<\J (f,t)\>=<\Jo\> (A9) 

with ph the number density of the carriers (holes in our case) . In Heisenberg 
representation, the operator equation of motion for J is identical in form to 
the classical one, i.e. 

^ = —JoE a + — V e abc J b B c (A10) 
at m a m a ^ 

b,c 

and ^j7"o = on account of the charge conservation. Taking the expectation 
value of (A. 10) with respect to | >, we find that 

dj a e 2 p h e ^ 

-7T = E a H > e abc j b B c . (All) 

at m a m a ^ 

b,c 

The anisotropy can be scaled away as: 



3a = (A12) 



'rrin 



E a = ^T a E a (A 13) 

and 

B c = ^m~m b B c (A 14) 

with a, 6, c a permutation of 1, 2, 3. Then 

^ = e 2 p h E a + e^2 e abc3bB c . (A 15) 

For E a = E a (u)e- iut , we find j a = j a {^)e' iu}t with 

j a (u>) = J2°ab{u)E b (u>), (A16) 
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where 



M") = z4^y (<U - ^B a B b ) - -LE*^ £ e a&c S c (A 17) 



I ^2 g2 £>2 

with cur = e\/ — -L + — 2 „ + ^, i, the cyclotron frequency. Undoing the 

13 V "T>2 m 3 m3mi m\m2 J m ^ o 

scale transformations (A. 12-14), we obtain the complex conductivity tensor, 
^abi^) = ^m a7nb ^ ^abi^] )• Substituting into the standard formula, we end up 
with the dielectric tensor with respect to the lattice frame 

£ab(u) = 5 a b . ; o 2T ( ^ab oB a B b + % \ e abc B c 

y/m a m b (uj 2 - u z B ) V u z ^m a m b uj ^ 

(A18) 



Appendix B 



In this appendix, we shall derive the improved Pauli term (3.32) by means 
of Foldy-Wouthuysen transformation of Dirac equation. The single electron 
Dirac equation in an ionic crystal and in an external radiation field reads 

= H*, (B.l) 

with 

H = -ipir ■ (V - ieA) + eA + V(r) ■ r + psm e , (5.2) 

where (Aq,A) is the four component gauge potential with Aq containing the 
self-consistent electrostatic potential of ion cores and electron clouds, V(r) is 
the exchange part of the self-consistent potential, ip(f) is a 4 x 1 spinor wave 
function and p\ and tj are the same 4x4 matrices introduced in the section 
2 (in terms of the standard notation, p3 = pyr = a). Equation (B.l) can 
be regarded a relativistic Hartree-Fock equation with a filled Dirac sea of 
negative energy levels and a Fermi sea of the positive levels. As we shall see, 
to the order we are interested, the details of Aq(t) and V(r) are unimportant 
provide that their magnitudes are nonrelativistic and they produce the same 
conducting band given by (2.28). 

For a nonrelativistic electron, the lower 2x1 spinor of ip is suppressed 
by a factor ~ v/c relative to the upper 2x1 spinor. Following Foldy and 
Wouthuysen [12], we decompose the Hamiltonian (B.2) into three part, 

H = p 3 m e + £ + £>, (fl.3) 
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where 

S = eA (f)+ P3 V(r)-T (BA) 
O = -ipir- (V -iel(f)). (B.5) 

Among sixteen basic 4x4 matrices, 1, pi, r^, piTj (i = 1,2, 3, j = 1,2,3), 
Pi? P2, P\T{ and p2^i couple to the upper 2x1 spinor and the lower one and 
will be referred to as odd operator, the rest of them does not couple the 
upper 2x1 spinor and the lower one and will be referred to as odd operator. 
Therefore p 3 m e and £ is even while O is odd. A systematic nonrelativistic 
approximation of a Dirac equation amounts to a set of successive unitary 
transformations 

= e" iS V (B.6) 

and 

H' = H- ie lS ^-e' iS , (B.7) 
dt v ' 

which eliminates the odd operators to the desired order of For the Dirac 

Tile 

Hamiltonian (B.2) and to the order the following three transformations 
serves the purpose: 

e -iS = e -iSs e -iS 2e -iS 1 

with 

Si = -i^-psO (B.9) 
2m e 

S 2 = ^(0-i[0,£}) (5.10) 



and 

/4 

P3 

"e 

The transformed Hamiltonian reads 



Ss = 8i p3 (3^ + ^ + i[£,0-i[0, £]] - i[6, £] - i[0, £]) . (B.ll) 



H' = p 3 m e + eA - — !— p 3 (V - ieA) 2 —p 3 r ■ B 

2m e m e 

ie ; ?.VxE + ^r.Ex(V-ieA)--^V.E 



Am 2 Am 2 &m 2 

;P3 {r • (V - ieA), {t ■ (V - iel),V(r) ■ f}}, (B.12) 



"e 

e 

8m 2 



where E = — ^ — VA\ and B = V x A. If we neglect the local magnetic field 
produced by the spin ordering, A is entirely due to the radiation field, so is B. 
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The electric field E, however, can be decomposed into a static field caused by 
the ionic potential and the radiation field. The Hamiltonian is a quadratic 
function of the radiation field, A, -Z? ra d. and B. For a nonrelativistic Fermi 
sea of plasma frequency u p and a photon of frequency u ~ u p , k ~ ui. With 
\V\ ~ uj p and eAo(r) ~ cu p , it is easy to see that the leading contributions 
to the antisymmetric part of the self energy function come from the second 
order perturbation of the Pauli term and the first order perturbation of the 
term 

~^2j-AxE iaA . (BA3) 
of (B.12). The latter gives rise to the Farady angle (3.28). 
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